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ABSTRACT 

We present an analysis of the intrinsic UV absorption in the Seyfert 1 galaxy Mrk 279 based on 
simultaneous long observations with the Hubble Space Telescope (41 ks) and the Far Ultraviolet Spec- 
troscopic Explorer (91 ks). To extract the line-of-sight covering factors and ionic column densities, we 
separately fit two groups of absorption lines: the Lyman series and the CNO lithium-like doublets. For 
the CNO doublets we assume that all three ions share the same covering factors. The fitting method 
applied here overcomes some limitations of the traditional method using individual doublet pairs; it 
allows for the treatment of more complex, physically realistic scenarios for the absorption-emission 
geometry and eliminates systematic errors that we show are introduced by spectral noise. We derive 
velocity-dependent solutions based on two models of geometrical covering - a single covering factor for 
all background emission sources, and separate covering factors for the continuum and emission lines. 
Although both models give good statistical fits to the observed absorption, we favor the model with 
two covering factors because: (a) the best-fit covering factors for both emission sources are similar for 
the independent Lyman series and CNO doublet fits; (b) the fits are consistent with full coverage of 
the continuum source and partial coverage of the emission lines by the absorbers, as expected from the 
relative sizes of the nuclear emission components; and (c) it provides a natural explanation for vari- 
ability in the Lya absorption detected in an earlier epoch. We also explore physical and geometrical 
constraints on the outflow from these results. 

Subject headings: galaxies: individual (Mrk 279) — galaxies: active — galaxies: Seyfert — ultraviolet: 
galaxies 



1. INTRODUCTION 

Mass outflow, seen as blueshifted absorption in UV 
and X-ray spectra, is an important component of 
active galactic nuclei (AGNs' see recent review in 
iCrenshaw. Kraemer. fc George I [20031) . This "intrinsic 
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absorption" is ubiquitous in nearby AGNs, appearing 
in over half of Seyfert 1 galaxies having high-quality 
UV sp ectra obtained with th e Hubble Space Telescope 
{HST; ICrenshaw et al.l 119991) and the Far Ultraviolet 
Spectroscopic Explorer (FUSE: Ikriss 1 1200^ . Spectra 
from the Advanced Satellite for Cosmology and Astro- 
physics (ASCA ) identified X-ray "warm absorbers" , seen 
as absorption edges, in a similar p ercentage of objects 
(jRevnolds ll997tlGeorge et al.ll998j) . Large total ejected 
masses have been inferred for these outflows, exceed- 
ing the accretion rate of the central black hole in some 
cases, indicating mass outflow plays an importan t role 
in the overall energetics in AGNs fe.g. iRevnolds l ll997'l. 
Recent studies have recognized and explored the poten- 
tial effect of outflows on all scales of the AGN environ- 
ment, from feeding the central supermassive black hole in 
AGNs (Blandford & Bcgelman _ 1999, 2004), to influcnc- 
ing the evolution of the h ost galaxy (Silk & Rees 199^ 
iScannapieco fc Oh I l2004 |) and the metallicity of the in- 
tergalactic medium l|Cavaliere. Laoi. fc Mem^l2002D . 

Measured ionic column densities provide the basis 
for interpretation of the physical nature of AGN out- 
flows. Detailed UV spectral studies over the past 
decade have shown measurements of these crucial pa- 
rameters are often not straightforward. Analyses of 
absorption doublets and multiplets has revealed the 
absorbers typically only partially occult the back- 
ground emission sources. Without proper treatment 
of this effect, the column densities could be severely 
in error fe.g. iWampler. Bergeron, fc Petitiean I 119931: 
iBarlow fc Sargent I Il997t iHamann et alJ Il997r Addi- 
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tional complications that could affect column density 
measurements are different coveri ng factors for differ- 
ent backgro und emission sources ijGangulv et alJ ITQQgt 
iGabel et al.l 1200 31. velocity-dependent covering factors 
fe.g. lAravetaLlllQflfll) . and inhomog eneous dis tributions 
of absorbing material Kool. Korista. fc Arav 2002.). 

Many recent investigations of AGN outflows have 
focused on intensive multiwavelength observations of 
Seyfert 1 galaxies. Seyferts are well suited for these stud- 
ies because they include the brightest AGNs in the UV 
and X-ray. The X-ray spectra contain the imprint of the 
bulk of the outflow's mass, which can now be deblended 
into individual absorption lines with the high-resolution 
capabilities of the Chandra X-ray Observatory (CXO) 
and XMM-Newton Space Observatory, allowing detailed 
study. The high quality UV spectra available with HST 
and FUSE provide a complimentary, precise probe of the 
complex absorption troughs. Due to the relatively nar- 
row absorption in Seyfert outflows, the important UV 
doublets and multiplets are typically unblended, allow- 
ing measurements of these key diagnostic lines. 

We have undertaken an intensive multiwavelength ob- 
serving campaign with HST/STIS, FUSE, and CXO to 
study the intrinsic absorption in the Seyfert 1 galaxy 
Mrk 279. Mrk 279 was selected for this study because 
of its UV and X-ray brightness and the rich absorp- 
tion spectrum in both bands, i ncluding unblende d and 
well- resolved UV doublets (see iScott et all 12004 here- 
after SK04) . Additionally, it has minimal contamination 
by Galactic absorption and a relatively weak contribu- 
tion from a narrow emission line region (NLR), both 
of which can complicate measurements of the absorp- 
tion properties. As part of a series of papers devoted 
to this campaign, we present here a detailed study of 
the UV absorption in the combined STIS and FUSE 
spectra. We develop a new approach for measuring the 
covering factors and column densities in the absorbers, 
making full use of the high quality far-UV spectrum. 
These measurements provide the foundation for subse- 
quent analysis and interpretation of the mass outflow in 
Mrk 279, and provide novel geometric constraints. In 
parall el papers, we present analysis of the X-ray spec- 
trum (|Costantini et al. ||2QQ4D. inhom ogeneous models of 
the UV absorption l|Arav et al.ll2004D . an d density diag- 
nosti cs based on Ov K-shell X-ray lines l|Kaastra et al.l 
I200I . In future papers, we will present photoioniza- 
tion models of the UV and X-ray absorption and anal- 
ysis of absorption variability. In the next section, we 
describe the HST and FUSE observations and present 
an overview of the absorption spectrum; in §3, we review 
the standard doublet technique for measuring intrinsic 
absorption and, together with an Appendix, discuss im- 
portant limitations of this method; the formalism of our 
fitting method and results for Mrk 279 are described in 
§4; in §5, the fits are interpreted, and implications for 
physical constraints on the outflow are explored; flnally, 
a summary is presented in §6. 

2. OBSERVATIONS AND THE INTRINSIC ABSORPTION 
SPECTRUM 

2.1. Simultaneous HST/STIS and FUSE Observations 
of Mrk 279 

The nucleus of Mrk 279 was observed for a total of 41 
ks (16 orbits) with the Space Telescope Imaging Spectro- 



graph (STIS) on board HST between 2003 May 13 - 18 
and for 91 ks with f(75'£' between 2003 May 12 - 14. The 
STIS observation used the E140M grating, which covers 
1150 - 1730 A, and was obtained through the 0".2 x 
0".2 aperture. The spectrum was processed with CAL- 
STIS v2.16, which removes the background light from 
each echelle order using t he scattered light model from 
iLindler fc Bowers"! 1)2000(1 . Low residual fluxes in the 
cores of saturated Galactic lines indicate accurate re- 
moval of scattered light: typical fluxes in the cores are 
± < 2.5% of the local unabsorbed continuum flux levels, 
and mean fluxes averaged over the absorption cores are 
< 3% of the noise in the troughs. The flnal spectrum 
was sampled in 0.012 - 0.017 A bins, thereby preserving 
the 6.5 km s^^ kinematic resolution of STIS/E140M. 

We found that the standard pipeline processing did 
not yield a fully calibrated spectrum due to two ef- 
fects: a) the echelle ripple structure, due to the char- 
acteristic efficiency of the d etector along each ord er, is 
not completely removed fsee lHeap fc Brown 1119971) . and 
b) the sensitivity of the MAMA detectors has degraded 
with time, and the change has not been incorporated in 
the pipeline for the echelle gratings. In order to cor- 
rect for these effects in the Mrk 279 spectrum, we used 
multiple spectra of the white dwarf spectrophotometric 
standard, BD-l-28 4122, one of which was taken close 
in time to our observation. First, a composite stellar 
spectrum of BD+28 41 22 composed of FO S and STIS 
data (Bohlin. Kickinso n. fc Galzetti ll2(K)ll) was used to 
flux-calibrate a 1997 STIS spectrum of BD-H28 4122 that 
does not exhibit the echelle ripple structure. This flux- 
calibrated spectrum was then used to correct a STIS 
spectrum of BD-l-28 4122, taken on 3 May 2003 with the 
same grating and aperture as our Mrk 279 observation, 
and which does show the same ripple structure seen in 
the Mrk 279 spectrum. These two corrections were per- 
formed by dividing the fiducial spectrum for each step 
by the comparison spectrum, fitting a polynomial to the 
result, and then multiplying the comparison spectrum 
by that polynomial. We used the polynomials from the 
second step described above to correct the spectrum of 
Mrk 279. To obtain a smooth correction for each or- 
der in the Mrk 279 spectrum, we applied averages of the 
polynomials corresponding to the four adjacent orders. 
Although we were able to remove most of the echelle rip- 
ple structure in this way, some lower amplitude residual 
curvature remains in some orders. However, the intrin- 
sic absorption features measured in this study are well- 
corrected. 

The FC/S'E; spectrum, obtained through the 30" x 30" 
aperture, covers 905 - 1187 A. The spectrum was pro- 
cessed with the current standard calibration pipeline, 
CALFUSE v2.2.3. The eight individual spectra obtained 
with FUSE, from the combination of four mirror /grating 
channels and two detectors, were coadded for all expo- 
sures. Mean residual fluxes measured in the cores of 
saturated Galactic lines are consistent with zero within 
the noise (i.e., standard deviation of the fluxes) in the 
troughs of these lines, indicating accurate background re- 
moval. The spectrum was resampled into ~ 0.02 - 0.03 
A bins to increase the signal-to-noise ratio (S/N) and 
preserve the full resolution of FUSE, which is nominally 
~ 20 km s-^ 

To place the FUSE and STIS spectra on the proper 
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TABLE 1 
Kinematics of Absorption 
Components in Mrk 279^^ 



Component Velocity FWHM 
km s"'^ km s"'^ 



1 85 40 

2 -265 50 
2a -290 30 
2b -325 30 
2c -355 65 

3 -390 20 

4 -460 20 
4a -490 65 

5 -550 30 



''Measurements for components 2, 
2a, 2c, and 4a are from the N V A1242 
line; components 1, 2b, 3, 4, and 5 are 
from Si III and C III. 

'^Radial velocity relative to the sys- 
temic redshift adopted for Mrk 279, 
z =0.305 



wavelength scale, we followed the procedure described in 
SK04. The centroids of prominent, unblended Galactic 
interstellar absorption lines were measured and used as 
fiducials in comparing to the Galactic 21 cm Hi line in 
the line-of-sight to Mrk 279 l|Wakker et alJl200H) . The 
lines measured in the STIS spectrum are consistent with 
the 21 cm Hi line within the measurement uncertainties 
and thus required no correction. The FUSE spectrum 
showed substantial shifts relative to the adopted stan- 
dard. Due to non-linear offsets in the wavelength scale, 
local shifts were measured and applied individually to 
each spectral region containing intrinsic absorption fea- 
tures. 

To normalize the absorption, we fit the total intrinsic 
(i.e., unabsorbed) AGN emission in the FUSE and STIS 
spectra over each intrinsic absorption feature. This was 
done empirically by fitting cubic splines to unabsorbed 
spectral regions adjacent to the features, at intervals of 
'--^ 5 A . We also derived models for the individual con- 
tributions of the different emission sources (continuum 
and emission lines), since they are required for our anal- 
ysis. For the continuum source, we fit a single power law 
(/a oc A~") to the observed flux at two widely separated 
wavelengths that are relatively uncontaminated with ab- 
sorption or line emission features, A = 955 and 1500 
A. After first correc ting for Galactic extinction (using 
the extinction law of iCardelli. Clavton. fc Mathis 1119891 
with E{B — V) — 0.016), we find a best-fit spectral index 
a = 1.6. This power law model matches the few other 
line- free regions of the UV spectrum, i.e. A « 1150 - 
1180, 1330, and 1390 A, to within a few percent. Thus, 
for the emission line model, we simply subtracted the 
continuum power law model from the empirical fit to the 
total emission. 

2.2. The Far- UV Absorption Spectrum 

The full far-UV spectrum from our FUSE and STIS 
observations is shown in Figure 1. The active nucleus 
in Mrk 279 was in a relatively high flux state during 
this epoch; the UV continuum flux was similar to a 
2000 January FUSE observation, and w 7 times stronger 



than in FUSE and STIS spectra obtained in 2002 May 
Full treatment of these earlier observations is given in 
SK04. Qualitatively, the intrinsic absorption spectrum 
in our new observations is similar to the earlier epochs 
(although some important variations were detected that 
wfll be the subject of a later study). Here, we give a 
brief overview, and refer the reader to SK04 for a more 
thorough phenomenological discussion of the absorption. 

Absorption from a range of ions is detected around the 
systemic velocity of the host galaxy, between v = —600 
to -1-150 km s"'^; we adopt the redshift for Mrk 279 from 
SK04, z = 0.0305 ± 0.0003. Normahzed absorption pro- 
files for some of the prominent lines are shown in Fig- 
ure 2. The absorption is seen to be resolved into mul- 
tiple distinct kinematic components at the resolution of 
STIS E140M and FUSE, revealing striking differences in 
the kinematic structure of different ions. Low- ionization 
species appear in several narrow components (see Si ill in 
Figure 2, but also Sin, Cii, Cm, and Niii in SK04 Fig- 
ures 7 - 14). However, the more highly ionized O vi, Nv, 
and C IV doublets, which are the primary UV signatures 
of intrinsic absorption in AGNs, are much broader and 
have different centroid velocities. The Lyman lines ex- 
hibit the kinematic structure of the low-ionization com- 
ponents, but also appear in the lower outflow velocity re- 
gion coinciding with the high-ionization lines, v w —300 
to -200 km s-^ 

We adopt the component numbering system from 
SK04, which was based on the kinematic structure in 
Ly/3. In Figure 2, dotted vertical lines mark the cen- 
troids of the components in SK04 that exhibit narrow 
absorption structure in low-ionization species in the cur- 
rent spectrum but which have no corresponding structure 
in the high-ionization CNO doublets. Centroids of the 
components seen in the high-ionization lines are identi- 
fied with dashed lines; we have added component 2c to 
the SK04 system based on structure in the C IV and N V 
profiles. Measured centroid radial velocities and widths 
of the components are listed in Table 1. The differences 
in ionization and kinematic structure between these two 
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Fig. 1.— Full far-UV spectrum of the active nucleus in Mrk 279 from the FUSE and HST/STIS spectra obtained in May 2003. Emission 
lines are labeled above the spectrum. The data were heavily binned for clarity in presentation. 



groups of components strongly suggests they are physi- 
cally distinct. Based on their ionization, narrow widths, 
distinct centroid velocities, and (in component 4) low 
density implied by the stringent upper limit on the C ii 
column density in the excited fine-structure level, SK04 
concluded at least some of the low- ionization components 
are not associated with the AGN outflow. Instead, they 
posited they are associated with gas located at relatively 
large distances from the nucleus - perhaps from an in- 
teraction with the companion galaxy MGG-l-12-13-024, 
high-velocity clouds associated with the host galaxy of 
Mrk 279, or, in the case of component 1, in the interstel- 
lar medium of the host galaxy. 

In this study, we restrict our attention to the bona fide 
intrinsic absorption, i.e., that presumed to be directly 
associated with an outflow from the AGN. We take this 
to include all absorption from the broad O vi, N v, and 
C IV doublets; Figure 2 shows that any absorption as- 
sociated with the narrow low-ionization components will 
at most only effect the outer wings of the outflow com- 
ponents in these lines. Conversely, Figure 2 shows Ly- 
man line absorption from the low-ionization components 
is strong and heavily blended with the intrinsic absorp- 
tion components; thus, we limit our analysis of H I to the 
uncontaminated region, v « —300 to —200 km s~^. 



3. THE DOUBLET METHOD: OVERVIEW AND 
LIMITATIONS 

We present here a brief review of the standard tech- 
nique for measuring UV absorption in AGN outflows and 
describe some limitations of this method to highlight the 
motivation for our method of analysis. In earlier stud- 
ies of intrinsic absorption, the red members of doublet 
pairs were often found to be deeper than expected rela- 
tive to the blue lines, based on their intrinsic 2:1 optical 
depth ratios. In many cases this was interpreted as due 
to partial coverage of the background nuclear emission by 
the ab s orbing gas, e.g.,IW amplcr. Bergeron. & Pctiticaiy 
l)1993|) . lBarlowfc Sargent. (,1997.'l . ,Hamann et al. l,1997j) : 
(other possibilities are scattering from an extended re- 
gion and emission from an extended source unrelated 
to the central engine of the AGN; ijC o hcn et al, 199^ 
iGoodrich fc Miller I [19951 iKraemer et alJl200lj) . If par- 
tial coverage is not accounted for, the absorption ionic 
column densities can be severely underestimated, which 
will dramatically affect the interpretation of the outflow. 
The expression for the observed absorption that includes 
the effects of line-of-sight covering factor (C) and optical 
depth (r) is: 



liv) ^ {1 - C{v)) + C{v)e-^^''\ 



(1) 
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Fig. 2. — Normalized absorption profiles from 2003 May STIS 
and FUSE spectra. The spectra are plotted as a function of 
radial velocity with respect to the systemic redshift of the host 
galaxy. The centroid velocities of kinematic components associ- 
ated with the AGN outflow are identified with dashed vertical 
lines. Components identifying low-ionization absorbers likely un- 
related to the outflow (see text) are shown with dotted lines. The 
difference in kinematic structure is evident in comparing the high- 
ionization CNO doublets with Sim. The hyS profile is not plotted 
at i; < —400 km due to contamination with Galactic absorp- 
tion at these velocities. 



where / is the normalized flux, and all quantities are 
written as a function of radial velocity, v. Since the opti- 
cal depths of the UV doublet pairs are in the simple 2:1 
ratio, equation 1 can be solved for the covering factor 
and optical de pths of each doublet IBarlow fc Sargentl 
[1997; Haman n et al.lll997fl . The resulting expressions, 
which we will refer to as the doublet solution, are: 



r. = -ln(^), (3) 

where r and b subscripts identify the red and blue mem- 
be rs o f the doublet, and the equation for r was derived 
bylArav, Korista, & dc Kool ( 2002). These expressions 
can be evaluated for unblended (i.e., sufflciently narrow) 
absorption doublets with members that are individually 
resolved, and derived as a function of radial velocity. 



While this has provided a revolutionary advance in the 
study of intrinsic AGN absorption, there are some key 
limitations to this method, described below. 

3.1. Multi- Component Nature of the Background 
Emission 

An implicit assumption in the doublet solution is that 
the absorption is imprinted on a uniform, homogeneous 
background emission source, since it allows for the solu- 
tion of only a single C and r. However, the AGN emission 
is comprised of multiple, physically distinct sources, i.e., 
a continuum source and emission line regions (includ- 
ing multiple kinematic components), which have differ- 
ent sizes, morphologies, and flux distributions. Thus, in 
cases where the absorber only partially occults the total 
background emission, the distinct sources would be ex- 
pected to have different line-of-sight covering factors, in 
general. This possibility was first explored by Ganguly et 
al. (1999) for the continuum source and broad emission 
line region (BLR), and was demonstrated i n the intrin- 
sic absorption systems by Ganguly et al., IGabel et alJ 
(2003), HaU ct al. (2003), and SK04. 

To account for multiple discrete background emission 
sources, equation 1 can be expanded to give the normal- 
ized flux of the j*'* line 

/, = E4i?}(qe--^ + i-q)], (4) 

where the i*'' individual emission source contributes a 
fraction i?J = /T,i[Fj] to the total intrinsic flux and 

has covering factor Cj. The effective covering factor for 
each line is the weighted combination of individual cov- 
ering factors, and can be written: 

Q = (5) 

These are expansions of the expressions given in 
iG_angulv et al. (1999) to include an arbitrary number of 
emission sources. 

The multi-component nature of the background emis- 
sion has several important implications for the analysis 
of AGN outflows: 

• From equation 5 it can be seen that lines of the 
same ion could have different effective covering factors, 
which may introduce an error into the doublet equation. 
This happens when the underlying emission fluxes dif- 
fer from the spectral position of one line to the other 
and is illustrated in the Nv doublet absorption shown 
in Figure 3. Here, the emission line flux underlying 
the blue member is « 15% greater than under the red 
line, while the continuum flux under the two lines is 
identical. The magnitude of this error depends on the 
slopes of the flux distributions between the doublet lines 
and differences in individ ual covering factors o f the dis- 
tinct emission sources. iGangulv et alJ l|1999r) showed 
this effect is typically small when considering the con- 
tinuum/BLR distinction, due to the gradual slope of the 
BLR; however, if the doublet lines are near saturation, 
it could have a very large effect since small flux differ- 
ences correspond to large optical depth differences in 
these cases. Also, an underlying narrow emission line 
component could have a pro nounce d effect on the so- 
lution (Arav, Koris^^^^^^ool 200 ^ [Kraemer et alJ 
2002: Gabel. Kraemer. fc Crensh aw 2004^! 

• Without separation of the covering factors of the dis- 
tinct sources, covering factors derived from a doublet pair 
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Fig. 3. — Spectrum of Lyo - Nv (top panel) and Ly/3 - Ovi 
(bottom panel), illustrating complexities in treating covering fac- 
tors in intrinsic absorption measurements. The continuum flux 
level, plotted as dashed lines, is essentially identical for each pair 
of lines, while the emission-line fluxes underlying each line differs 
greatly. This has important implications for the relative effective 
covering factors. Additionally, the nature of the BLR emission 
underlying these lines is complex: Ly/3 absorption lies on the high- 
velocity blue wing of O VI, while Nv absorption has a contribution 
from the red wing of the Lya BLR. 



cannot be applied to measure column densities of other 
lines. This is evident in Figure 3; clearly the effective cov- 
ering factor derived for N v is not appUcable to Lya if the 
individual continuum and emission-Hne covering factors 
differ since Lya has much more underlying line flux. Sim- 
ilarly, Ly/3, plotted in the bottom panel in Figure 3, will 
not generally have the same effective covering factor as 
Lya, due to the different emission line contributions un- 
der each line. However, if the individual covering factors 
are known, effective covering factors can be constructed 
for any line using equation 5. This is important for mea- 
suring singlet lines or contaminated multiplcts that have 
no independent measure of the covering factor. 

• Finally, the doublet solution misses potentially im- 
portant, unique constraints on the absorption and emis- 
sion geometry. For example, combined with estimates 
of the sizes of the individual sources derived from other 
techniques, the individu al covering factors constrain the 
size of the absorber, e.g. iGabel et al.l 1)2003(1 . and the rel- 



ative location of the different emission components and 
absorber as projected on the plane of the sky. The indi- 
vidual covering factors can also serve as a unique probe 
of more detailed geometry of the background emission. 
Consider for example the O vi - Ly/3 spectrum shown 
in Figure 3. The Ovi doublet absorbs its own emis- 
sion line flux at the blueshifted velocity of the outflow 
{vblr ~ —600 to — 200 km s~-^), while the Ly/? absorber 
sits primarily on the high- velocity blue wing of the O vi 
BLR at vblr ~ —2200 km s~^. The Ly/3 line emission 
is relatively weak. Similarly, there is a contribution from 
the extreme red wing of the Lya BLR profile under the 
Nv absorption. Constraints on the emission line cov- 
ering factors for these lines could be used to probe the 
kinematic-geometric structure of the BLR; the absorber 
can thus serve as a filter to view and explore the back- 
ground AGN sources. 

3.2. Systematic Errors in the C ~ t Solutions 

Another limitation is that the doublet method always 
gives a solution, but it is often difficult to gauge its ac- 
curacy due to the non-linear dependency of the solution 
on measurement errors. To explore this, we have gen- 
erated synthetic absorption profiles that include random 
fluctuations simulating spectral noise, and calculated C 
and T using the doublet equations. Illustrative results 
are shown in an Appendix, whereas a complete quanti- 
tative treatment will be presented in a later paper. We 
find there are systematic errors in the solutions that can 
give misleading results. These errors are not random 
about the actual value, but rather systematically under- 
estimate the actual covering factor, as seen in the Ap- 
pendix; the discrepancy in the solution increases with 
weaker absorption doublets a nd decreased S/N. Addi- 
tionally, 'GaneiHT'et'ar (1999) demonstrated that the fi- 
nite instrumental line spread function can lead to further 
systematic errors in the doublet solution. 

4. OPTIMIZATION FITTING OF THE INTRINSIC 
ABSORPTION: LYMAN SERIES AND CNO DOUBLET 
GLOBAL LINE FITS 

Motivated by the limitations of the traditional doublet 
method described above, we introduce here a different 
approach for measuring intrinsic absorption. The under- 
lying principle is to increase the number of lines that 
are simultaneously fit in order to (a) explore additional 
parameters contributing to the formation of observed ab- 
sorption troughs and (b) overconstrain the set of equa- 
tions. This allows for the treatment of more complex, 
physically realistic scenarios of the absorption-emission 
geometry. By minimizing errors to simultaneous fits of 
multiple lines, noise in the spectrum will generally be 
smoothed out, in contrast to the erratic behavior of the 
doublet solution demonstrated in §3.2 and the Appendix. 

4.1. Formalism 

Our fitting algorithm employs the Levenberg- 
Marquardt non-linear least-squares minimization 
technique to solve equation 4 for specified absorption 
parameters (Cj, Tj)^'^. It is similar in principle to 

Using software provided by C. Markwardt, 
http://cow.phys.wisc.edu/~craigm/idl/idl.html] which is based 
on the MiNPAGK-1 optimization software ot J. More available at 
www.netlib.org 
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that used in SK04 to analyze the Lyman lines in 
earlier spectra of Mrk 279. Given a total of n observed 
absorption lines (Ij) as constraints, up to n— 1 param- 
eters can be modeled. No a priori assumptions are 
made about the kinematic distribution of the covering 
factors and optical depths of the absorbing material 
(e.g., Gaussian) - indeed, one goal is to solve for the 
velocity-dependent absorption parameters to constrain 
the kinematic-geometric structure of the mass outflow. 
Thus, we derive fits to the absorption equations for each 
velocity bin. This also avoids errors in the solutions 
resulting from averaging over variable profiles. The 
algorithm minimizes the function, with each data 
point appropriately weighted by the 1 a errors, which 
are a combination of spectral noise and estimated 
uncertainties in fitting the intrinsic underlying fluxes. 
For the latter, the continuum flux uncertainties were 
determined from the residuals between the power-law 
model and the line-free regions of the spectrum. We 
estimated uncertainties in the emission-line fluxes by 
testing different empirical fits over the absorption 
features, finding the range that gave what we deemed 
reasonable line profile shapes. 

The key requirement in employing this technique is to 
link multiple absorption lines for simultaneous fitting. 
There are two general ways to do this: 

• Lines from the Same Energy Level: The most 
straightforward way is to fit all available lines arising 
from the same ionic energy level, thereby eliminating 
uncertainties in ionic abundances or level populations. 
If the relative underlying fiuxes from distinct emission 
sources differs between the lines, the individual cover- 
ing factors of those emission sources can be derived. In 
the subsequent analysis, we fit all of the uncontaminated 
Lyman series lines in Mrk 279. These lines are ideal 
because they span a very large range in optical depth 
and have significantly different amount s of underlying 
emission- line flux fsee iGabel et aLll2003j) . Additionally, 
the full set of lines is accessible in low redshift AGNs with 
combined FUSE and STIS spectra. We note the Fell UV 
multiplets,_which appear in a small fraction of AGN ab- 
sorbers (|de Kool et al.l 120011 iKraemer et aI]l200H) . are 
another promising set of lines for this analysis. 

• Global Fitting Approach: The second approach in- 
volves linking lines from different ions (or any group of 
lines arising from different levels) by placing physically 
motivated constraints on their absorption parameters. In 
our analysis of Mrk 279 below, we fit the six combined 
lines of the Ovi, Nv, and CiV doublets by assuming 
they share the same covering factors. Another potential 
application of this method is to link the absorption for 
a given line in spectra from different epochs via assump- 
tions about the relative values of the absorption parame- 
ters between epochs (e.g., assuming the covering factors 
did not change). The validity of the assumptions used 
to link the equations can then be tested by the result of 
the fit. 

4.2. Covering Factor and Optical Depth Solutions for 

Mrk 279 

For the intrinsic absorption in Mrk 279, we indepen- 
dently fit the two groups of lines described above: the 
Lyman series lines and the combined CNO doublets, i.e. 
global fit. For each set of lines, we tested two differ- 



ent models of the absorption covering factor. In model 

A, a single covering factor was assumed to describe all 
lines, i.e., no distinction was made between the different 
emission sources. In model B, independent covering fac- 
tors for the continuum source and emission lines, and 
C', were assumed. In this case, the gene ral expressions 
in equ ations 4 and 5 reduce to those in iGangulv et all 
llT999l) . 

For the Lyman lines, the solvable range is limited to 
—300 ^ V < —200 km s^^, due to blending with the 
narrow, low-ionization components in the high velocity 
region of the outfiow (see §2.2). This is due to the failure 
of equation 4 where multiple absorption components with 
different covering factors contribute in the same velocity 
bin; there is no straightforward way to disentangle how 
the different absorbers overlap as projected against the 
background emission sources. Figure 2 shows Lya, Ly/3, 
and Ly7 give the best constraints for the Lyman series 
analysis, exhibiting clean well-defined absorption profiles 
at relatively high S/N. Weak absorption in LyS is also 
present in component 2, while Lye is contaminated with 
Galactic H2 absorption and thus omitted from the fitting. 
Ly^ is not detected in components 2 - 2a within the 
limits of the spectral noise, thus to reduce the effect of 
noise on the solution, we set the normalized fiux to unity 
for this line. All lines of higher order than LyC were 
omitted from the analysis since they exhibit no intrinsic 
absorption and provide no additional constraints. Thus, 
there are five lines as constraints to fit the two and three 
free parameters for the Lyman solution in model A and 

B, respectively. 

For the glolDal fitting of the Ovi, Nv, and CiV dou- 
blets, the absorption equations were linked by assum- 
ing their covering factors are equal (separately for the 
continuum and emission lines in model B). The optical 
depth of each ion is a free parameter. As described in 
§2.2 and seen in Figure 2, these lines are not strongly 
contaminated by the narrow, low-ionization absorption 
that dominates the Lyman lines; at most, there is only 
weak contamination in the wings of the broad, intrinsic 
absorption features by these low-ionization systems. 

Best-fits to the covering factors and optical depths (or 
equivalently column densities) for both line groups are 
shown in Figure 4 (model A) and Figure 5 (model B). 
In all cases, the parameter space search for the covering 
factors was restricted to the physically meaningful range, 

< C* < 1. The plotted error bars represent the formal 

1 a statistical errors in the best-fit parameters. Specif- 
ically, these correspond to values giving Ax^ =1, and 
were computed from the diagonal elements of the covari- 
ance matrix for the optimal fit (Bcvington 1969). For 
cases where the covering factor solution is a boundary 
value (0,1), e.g. much of the C"^ solution for the CNO 
doublets in model B, no covariance matrix elements are 
computed for that parameter since it is not a minimum 
in the solution. For these cases, we estimated uncertain- 
ties by deriving solutions for models keeping the param- 
eter fixed, and finding the value giving A^^ =1 from the 
best-fit solution at the boundary value. To ensure the 
computations did not erroneously stop at local minima, 
we generated solutions using different starting points for 
the parameter search space; identical results were found 
in all cases. The fitted profiles, derived by inserting the 
best-fit solutions into equation 4, for both models (model 
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Fig. 4. — Best- fit covering factor and column density/optical 
depth solutions from minimization for the single covering factor 
geometric model (model A). Top panels show C solutions, solved 
independently for the Lyman series (left panels) and combined 
CNO doublets (right). Bottom panels show the Hi column density 
and Civ (dashed histogram), Nv (dotted), and Ovi (solid) optical 
depth solutions. The contaminated, and thus unreliable, region of 
H I absorption is plotted with dotted histograms. 
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A in blue, model B in red) are compared with the ob- 
served normalized absorption profiles in Figure 6. 

4.3. Errors and Uncertainties 

There are some possible errors in this fitting that 
should be mentioned. First, the emission lines are treated 
as all arising from a single component. Line emission 
from distinct kinematic components, i.e. broad, interme- 
diate, narrow line regions, that are covered by different 
fractions by the absorbers could introduce an error into 
the solution (although, the narrow line region in Mrk 279 
is relatively weak). Also, there are cases where absorp- 
tion features sit on the BLR emission from different lines 
and thus sample different velocities of the BLR profile 
(i.e., Ly/3 and Nv described in §3.1 and shown in Figure 
3). This could introduce an error into the solution if there 
are spatial inhomogencitics in the BLR g function 
of velocity. Second, the area on the sky sampled by the 
FUSE aperture is four orders of magnitude greater than 
the STIS aperture. While this has no consequences for 



the continuum and BLR emission, which are unresolved 
and much smaller than the 0."2 STIS aperture, any ex- 
tended emission might effect the solutions. This could in- 
clude scatt ering of nuclear emiss ion by an extended scat- 
terer ('e.g. iKraemer et all 120011) or extended Ovi NLR 
emission. Finally, we have assumed the absorption opti- 
cal depth of each line (r^) is uniform across the lateral 
extent of the emission source s, and thus the absorber is 
completely homogeneous fsee lde Kool. Korista. &: Aravl 
l2002tl . Models departing from this assumption are pre- 
sented in lArav et al. (20041 . 

5. DISCUSSION AND INTERPRETATION 

5.1. Favored Absorption Model: Independent 
Continuum and Emission Line Covering Factors 

Figure 6 shows both geometric models are able to 
match the intrinsic absorption features well at most out- 
flow velocities, indicating the solutions for these models 
and parameters are degenerate over much of the profiles. 
There are some regions fit somewhat better by model B, 
particularly CiV in the low- velocity outfiow component. 
However, there is additional, stronger evidence that sup- 
ports the two covering factor model for the outfiow, as 
outlined below. 

• Consistent Covering Factor Solutions From the In- 
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Fig. 5. — Same as Figure 4 for the two covering factor geometric 
model (model B). Continuum covering factors are plotted in black 
and emission-line covering factors in red. 
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dependent Lyman Line and Global CNO Doublet Fitting 
Methods: Figure 5 shows that both the Lyman and global 
CNO doublet best-fit solutions are consistent with full 
coverage of the continuum source over most of the veloc- 
ity range with a valid Lyman solution. Further, the inde- 
pendent solutions for the two groups of lines are plotted 
together in Figure 7 for a direct comparison; the emis- 
sion line covering factor profiles computed in model B are 
shown in the top panel and the single covering factors de- 
rived in model A are in the bottom panel. The H i and 
CNO C' fits are nearly identical over most of the core 
region of the absorption, with « 0.7 (although in the 
wing of the absorption, v > —250 km s~^, the solutions 
diverge somewhat). In contrast, the Lyman method so- 
lution in model A is systematically less than the global 
doublet solution by about 0.1 - 0.15 at all velocities. 

• Consistency with the Emission Source Sizes: The 
continuum and emission line covering factor solutions in 
model B are also physically consistent with our under- 
standing of the sizes of those emission sources. Reverber- 
ation mapping studies of AGNs have shown the BLR is 
substantially larger than the ionizing continuum source. 
These studies measure the BLR to be several to tens of 
hght-days across (e.g. iPeterson et aLllTOQ^ iKaspi et al.l 



|2000() . while the UV continuum source is likely at least 
an order of magnitude small er in size and possibly much 
smaller fe.g. iLaor fc Netzer I 1989: Proga 2003). Thus, 
the solution of a fully covered continuum source and par- 
tial coverage of the emission lines is consistent with the 
nuclear emission geometry. The fact that this result was 
arrived at separately with the independent Lyman and 
global CNO fits is compelling. 

• Variability in Lya Absorption: In Figure 8, the nor- 
malized Lya profile from our observation is compared 
with the STIS spectrum obtained a year earlier, showing 
the absorption was shallower in the previous epoch by ~ 
0.15 in normalized flux units. This variability is not due 
to a lower H I column density; the strength of Ly/3 in a 
contemporaneous FUSE observation indicates Lya was 
saturated in the earlier epoch (SK04), as it is in the 2003 
spectrum, hence the profiles in both epochs simply de- 
lineate the unoccultcd flux. Thus, a change in covering 
factor must be responsible for the observed variability. 
Noting that the emission line-to-continuum flux ratio was 
greater in the 2002 spectrum than in 2003, we tested if 
the difference in covering factor of the respective sources 
could explain the variability. To this end, we constructed 
a synthetic absorption profile for Lya in the 2002 epoch 
based on the results from our analysis of the 2003 spec- 
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Fig. 6. — Best-fit intrinsic absorption profiles for tlie Lyman lines 
(left) and CNO doublets (right) from the two geometric models. 
The normalized observed spectrum is plotted in black, model A 
profiles are in blue, and model B profiles in red. The model profiles 
were derived from the best-fit covering factors and column densities 
shown in Figures 3 and 4 using equation 4. The contaminated 
region of the Lyman line fit (dotted lines in Figures 3,4) is not 
plotted. 



trum. Specifically, we set C"^ =1 at all velocities, and 
solved for C' in the 2003 Lya profile. The observed con- 
tinuum and emission line fluxes in the 2002 spectrum 
were then weighted by these covering factors according 
to equation 4. The resulting model profile, shown in the 
bottom panel of Figure 8, matches the observed absorp- 
tion very well. This provides a natural explanation for 
the variability in the Lya absorption profile: the cover- 
ing factors of each source individually remain the same, 
but a change in the effective covering factor occurs due 
to different relative strengths of the distinct background 
emission components. 

5.2. Comparison of Solutions from Different Methods 

The covering factor and optical depth solutions for 
Civ, Nv, and Ovi from the two-C global fit (model B) 
are compared with the single doublet solutions (equa- 
tions 2 and 3) in Figure 9. The global fit effective cover- 
ing factors are weighted combinations of the individual 
covering factors shown in Figure 5, derived using equa- 
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tion 5. In Figure 9, solid red lines show C profiles for the 
short-wavelength members of each doublet and dotted 
lines show results for the long-wavelength lines; they are 
nearly identical in all cases, indicating negligible errors in 
the doublet solution due to different contributions from 
underlying line emission (see first bullet in §3.1). 

Figure 9 shows the solutions from the two methods dif- 
fer significantly in some regions, with important implica- 
tions for the kinematic-geometric and -ionization struc- 
ture in the outflow. For example, the doublet cover- 
ing factor solution, plotted in black, exhibits much more 
velocity-dependence than the global fit. This is seen, for 
example, in the wings of Civ and Nv in component 4a 
and the red wing of component 2, where C decreases 
strongly in the wings, while the optical depth profiles 
show little or no relation to the absorption trough struc- 
ture when compared with Figure 2. Further, the kine- 
matic structure in these ions is much different from the 
O VI solution. In these cases, the doublet solution implies 
the observed absorption profiles are determined almost 
entirely by velocity-dependent covering factor. Addition- 
ally, in some regions, particularly in C IV, the doublet so- 
lutions for C and t exhibit sharp fluctuations over small 
velocity intervals that do not coincide with structure in 
the observed absorption features, nor with the solutions 
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from the other ions (e.g., components 2a - 2c, and 4a). 

Comparison with the observed profiles in Figure 2 
shows these effects in the doublet solution are strongly 
correlated with the absorption strength. Smaller derived 
covering factors are associated with weaker absorption, 
both in the wings of individual kinematic components 
and in overall absorption features, such as C iv compo- 
nents 2c and 2a. As discussed in the Appendix and §3.2, 
synthetic absorption profiles show the doublet solution 
suffers from systematic errors consistent with this general 
trend - spectral noise causes an underestimate of the cov- 
ering factor. The simulations we present in the Appendix 
indicate the solutions become unreliable for cases where 
T < 1 in the blue doublet member, and increasingly so for 
lower r. This, combined with the large, seemingly ran- 
dom fluctuations seen in certain regions suggests some 
features in the doublet solutions in Figure 9 are artifacts 
due to noise and not real features in the outflow. How- 
ever, one puzzling aspect of the solutions, in comparison 
to the results in the Appendix, is that there are not more 
negative values for the covering factors. The simulations 
predict that for sufficiently low optical depth, there is a 
high probability that C <0. This should be especially 
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Fig. 7. — Comparison of covering factor solutions for the two in- 
dependent methods — Lyman and global CNO doublet lines. Best- 
fit solutions to the emission-line covering factors from model B are 
shown in the top panel. The single covering factor solutions from 
model A are in the bottom panel. 
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Fig. 8. — Variability in Lya absorption. The top panel shows Lya 
absorption in components 2-2a was shallower in a May 2002 STIS 
spectrum (solid histogram) than in May 2003 (dashed). The bot- 
tom panel shows a model of the 2002 profile (dotted), constructed 
with individual covering factors derived from the 2003 spectrum 
(Figure 5), matches the observed profile well. This indicates the 
variation in absorption depth is likely due to a change in effective 
covering factor resulting from different relative strengths of the 
continuum and emission-line fluxes in the two epochs. 



apparent for intermediate values of r, such as the middle 
two models presented in the Appendix, where dramatic 
fluctuations between large positive and negative values 
are expected. It is unclear why this is not more preva- 
lent in the solutions in Figure 9. 

In Table 2, we compare integrated column densities 
measured from the different solutions: the two geometric 
coverage models presented in §4.2, plus the traditional 
single doublet solution for the CNO doublets. The in- 
tegrated column densities for the doublet solution are 
larger than the global fits in all cases. In component 
2+2a, the doublet solution is 60% greater than the 2-C 
model for C iv, while in component 4a, it is 45% and 
70%i greater for Civ and Nv, respectively. Comparing 
the two global fit models, the CNO column densities from 
model B are greater than the single C model by - 30% 
, while the Hi column density is 50% less. 
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Fig. 9. — Comparison of covering factor and optical depth solu- 
tions derived with the global-fit (red) and single doublet (black) 
methods. Left panels: Covering factors from the global-fitting 
method are the effective values derived by weighting the individ- 
ual covering factors in model B by their respective fluxes. The 
short-wavelength doublet members are plotted with solid red lines 
and the long-wavelength members with dotted lines. Right panels: 
Optical depths for the long-wavelength lines of the doublets. 



5.3. Constraints on the Nuclear Absorption - Emission 

Geometry 

Here we explore constraints on the nuclear absorption 
and emission geometry available from the covering fac- 
tor solutions derived in §4.2. Full analysis of the physical 
conditions in the outflow, using photoionization model- 
ing of the combined UV and X-ray absorption, will be 
presented in a future study. 

As discussed in §5.1, we adopt the 2-C global fit 
(model B) for the covering factors and column densi- 
ties. The covering factor solutions in Figure 5 are con- 
sistent with the UV continuum source being fully in our 
sight-line through the absorber, for all kinematic compo- 
nents, while the emission lines are only partially so. This 
constrains the relative line-of-sight geometry of the ab- 
sorbers and nuclear emission sources and places a lower 
limit on the transverse size of the outflow. Monitoring 
of time lag s between B LR and continuum variations in 
Mrk 279 bv lMaoz et all (499Q.1 gives the size of the BLR, 
Rblr = 12 Hght days. Thus, the projected size of the 
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UV absorbers on the plane of the sky, against the AGN 
emission, is at least (C')^/^ x Rblr, or «10 hght days 
in the cores of the absorption components, and possibly 
decreasing in the wings. 

An additional probe of the nuclear geometry is possible 
because the Ly/? absorption lies on the high- velocity blue 
wing of the O vi BLR profile, at vblr ~ -2200 km s~^, 
as discussed in §3.1 and illustrated in Figure 3. Thus, 
it samples different BLR kinematics than most other ab- 
sorption lines, which absorb BLR velocities coinciding 
with the absorption outflow velocities {vblr ~ —300 
km s~^). Therefore, comparison of the Ly/3 emission- 
line covering factor and those associated with other lines 
serves as a probe of the kinematic-spatial structure of the 
BLR. The potential affect of complex velocity-dependent 
structure in the BLR on the absorption covering factors 
was explored in Srianand & Shankaranaravanan ( 1999). 
Since we assumed a priori that all lines share the same 
individual covering factors in the Lyman series fit in §4.2, 
a legitimate comparison requires obtaining independent 
constraints on the covering factor of the O vi BLR blue 
wing by the Ly/3 absorber. This comes straightforwardly 
from the observed Ly7 absorption line, since Ly7 has no 
underlying line emission. Using the result that C"^ =1, 
the intrinsic optical depth ratio for Ly/5 : Ly7, and the 
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TABLE 2 

Integrated Ionic Column Densities'' 



Model A: 1 C Fit Model B: 2 C Fit Doublet Method Full Coverage'' 
Ion 2+23"= Aa" 2+2a 4a 2+2a 4a 2+2a 4a 



C IV 1.6 0.9 2.1 1.1 3.3 l.e'' 1.5 0.9 

N V >6.6 2.1 >6.6 2.1 >8.1 3.6 4.7 2.0 

O VI >14 >13 >14 >13 >14 >13 9.8 9.5 

HI 7.5 • • • 5.0 • • • • • • • • • 4.4« 



"AH column densities in units of 10^'' cm for doublets, measurements are for the 
red member; lower limits are quoted for lines with regions having r > 3 
=1 assumed everywhere. 
•^Component 2+2a, v = —230 315 km s~^; Component 4a, v = —435 535 

km s~^ 

Integrated over limited range due to infinities in solution. 
''Using Ly7 line. 



observed Ly7 profile, the absorption profile for Ly/j can 
be derived as a funetion of C' from equation 4. Some 
illustrative results are shown in Figure 10. This shows, 
for example, that both an unocculted and fully occulted 
high velocity Ovi BLR are ruled out by the data. We 
find that values ranging from 0.5 ^ C < 0.8 are required 
to fit the majority of the Ly/3 absorption in components 
2 - 2a. This is similar to the emission line covering factor 
derived from the global fit to the CNO doublets (Figure 
5b). Additionally, the Lya C' profile can be derived inde- 
pendently for comparison with Ly/3, in a similar manner 
as for Ly/3 (i.e., using C° =1, the observed Ly7 profile, 
and the Lya : Ly7 r ratio). The result is identical to the 
solution to the combined Lyman lines in components 2 
- 2a shown in Figure 5a within uncertainties. Therefore, 
the absorption covering factor of the Ovi BLR emis- 
sion at vblr ~ —2200 km s^^ is similar to C' at lower 
emission-line velocities (by the CNO doublets and Lya). 
These results may provide constraints for models of the 
BLR, e.g., testing disk vs spherical geometries and out- 
flow vs rotational kinematics for the BLR. 

6. SUMMARY 

We have presented a study of the intrinsic UV absorp- 
tion in the Seyfert 1 galaxy Mrk 279 from an analysis of 
combined long observations with HST/STIS and FUSE. 
These spectra were obtained simultaneously in May 2003 
as part of an intensive multiwavelength observing cam- 
paign. 

We present a review of the standard technique for mea- 
suring intrinsic UV absorption parameters based on indi- 
vidual doublet pairs, showing some key limitations of this 
method: 1) It cannot treat multiple background emis- 
sion sources. This introduces a potential error in the 
solution and misses important geometric constraints on 
the outflow. 2) Using synthetic absorption profiles, we 
show it systematically underestimates the covering fac- 
tor (and overestimates r) in response to spectral noise. 
The discrepancy in the solution is shown to be strongly 
dependent on absorption strength. 

To measure the UV absorption parameters in Mrk 279, 
we independently fit two groups of lines: the Lyman se- 
ries lines and the combined CNO lithium-like doublets. 
The doublet fitting involved a global fitting approach, 
which assumes the same covering factors apply to all 



ions. By increasing the iniinber of lines that are simulta- 
neously fit, more complex and physically realistic models 
of the absorption-emission geometry can be explored. So- 
lutions for two different geometrical models, one assum- 
ing a single covering factor for all background emission 
and the other separate covering factors for the continuum 
and emission lines, both give good statistical fits to the 
observed absorption. However, several lines of evidence 
support the model with two covering factors: 1) the inde- 
pendently fit Lyman lines and CNO doublets give similar 
solutions to the covering factors of both emission sources; 
2) the fits are consistent with absorbers that fully occult 
the continuum source and partially cover the emission 
lines, consistent with the relative sizes of the emission 
sources; and 3) observed variability in the Lya absorp- 
tion depth can be explained naturally by this model as a 
change in effective covering factor resulting from a change 
in the relative strengths of the emission components. 

Comparison of the traditional solutions based on in- 
dividual doublets and the global-fit solutions shows the 
former exhibits much stronger velocity dependence. This 
is seen as decreases in covering factor in the wings of in- 
dividual kinematic components, and as peculiar fluctua- 
tions in both r and C in other regions of relatively weak 
absorption. In light of the systematic errors shown to be 
inherent in the individual doublet solution, we conclude 
some of these effects arc likely artifacts of the solution 
and should be interpreted with caution. 

The covering factor solutions from our global fit con- 
strains the relative line-of-sight geometry of the ab- 
sorbers and nuclear emission sources. The derived emis- 
sion line covering factor, combined with the size of the 
BLR, constrains the projected size of the absorber to be 
> 10 light days. We utilize the coverage of the high veloc- 
ity O VI BLR by the Ly/3 absorber to explore kinematic 
structure in the BLR; we find no evidence for depen- 
dence of the absorber's BLR covering factor on the BLR 
velocity. 

Support for this work was provided by NASA through 
grants HST-AR-9536, HST-GO-9688, and NAGS 12867 
and through Chandra grant 04700532. We thank D. 
Lindler and J. Valenti for their assistance in correcting 
the STIS spectrum and C. Markwardt for making his 
software publicly available. We also thank the referee P. 
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Fig. 10. — Constraints on the covering factor of the high blueshifted velocity Ovi BLR by the Ly/3 absorber. The observed normahzed 
Ly/3 profile (black) is compared with models derived for three values of emission-line covering factor (red). Models were derived using the 
observed hyy absorption, as described in the text. Models with an unocculted (dashed) and fully occulted (dotted) BLR are ruled out. An 
emission line covering factor of 0.7 (solid), consistent with the global fits, matches most of the observed profile well. 
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study. 

7. APPENDIX 

Here, we address how the expressions for covering fac- 
tor and optical depth from the doublet solution (equa- 
tions 2 and 3) depend on noise in the spectrum. We 
have generated synthetic absorption profiles for doublet 
pairs that include random fluctuations to simulate spec- 
tral noise. We derive the C - r solutions for these syn- 
thetic profiles to determine any trends in the solution as 
a function of noise level and absorption strength. 

Our synthetic profiles were derived in the following 
way. The optical depth profiles were assumed to be Gaus- 
sian, parameterized by the width (tr) and peak optical 
depth in the core of the blue line {Tmax)- The "true" 
normalized absorption profiles for the lines are then de- 
rived from equation 1, with the covering factors set at 
a constant value across the profile for each doublet pair, 
and T{v)r = T{v)b/2 at all velocities. Thus, we have as- 
sumed in effect a single background emission source to 
avoid the complications described in §3.1. For compari- 
son with our study of Mrk 279, we have set the velocity 
resolution of our synthetic profiles to that of the STIS 



E140M grating, and the absorption width (cr = 50 km 
s~^) to be approximately consistent with intrinsic ab- 
sorption components 2-2c. To simulate spectral noise, 
we generated a normally distributed random number as- 
sociated with each velocity bin in the synthetic spectra. 
The noise level was normalized by selecting the desired 
S/N in the unabsorbed continuum and then weighting 
the noise by flux level, S/N(?;) oc /(w)-^/^, according to 
Poisson statistics. 

We generated profiles for a range of S/N, C (real), and 
Tmax- Here, we give some brief illustrative results, while 
reserving a full analysis for a later study. The left pan- 
els in Figure 11 show synthetic profiles for doublet pairs 
with Tmax = 2, 1, 0.5, and 0.25. In all cases, C —0.8 
and the noise level was normalized to be 3% in the con- 
tinuum. The middle two panels show the corresponding 
covering factor and optical depth profiles derived directly 
from equations 2 and 3, with the actual values marked 
with dashed lines for comparison. Errors in the doublet 
solution are immediately apparent. The covering factors 
are systematically underestimated, and the magnitude 
of error is strongly dependent on absorption strength. 
This is seen both in the lower covering factors derived in 
the cores of features with lower Tmax (means and stan- 
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dard deviations measured over the central 150 km s~^ are 
printed in each plot), and in the decrease in C computed 
in the wings of each profile. 

These systematic errors are due to non-linear effects in 
the doublet equations. This is seen most clearly by com- 
paring the numerator, = (/,. — 1)^, and denominator, 
D = lb — 2Ir + 1, in the expression for covering factor 
(equation 2). The right panels of Figure 11 show the val- 
ues of D and N/C (in red), which would be identical in 
each velocity bin for infinite S/N. Due to the forms of N 
and D, these quantities have very different dependences 
on noise; AN/N becomes much smaller than AD/D for 
weak absorption and, at sufBciently small r, N is less 
than the noise level of D. As decreases relative to the 
noise in D for weaker absorption, the probability that 
< D < N/C becomes vanishingly small, and the aver- 
age value of N/ D becomes increasingly small. 

These errors could have pronounced effects on the in- 
terpretation of the outflow. Each solution that underes- 
timates C overestimates r. Thus, ionic column densities 
are systematically overestimated, with increasing relative 
discrepancies in weaker doublets, leading to errors in de- 
termining the ionization structure and total gas in the 
absorber via photoionization models. Additionally, the 
errors in covering factor solutions will effect geometric 
inferences. For example, due to the high- ionization state 
of AGN outflows, lower-ionization species appearing in 
UV spectra are generally weaker. Thus, the increasing 
discrepancy in weaker absorption doublets may lead to 
the misinterpretation of ionic-dependent covering factors. 
Also, weaker absorption in the wings of an absorption 
feature could lead to apparent velocity-dependent cover- 
ing factors that are instead due to optical depth varia- 
tions, or at least exaggerate the effect. 
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Fig. 11. — Systematic errors in the doublet solution resulting from spectral noise. Left panels show synthetic normalized doublet profiles 
(long wavelength members in red) with simulated noise for a range of absorption strengths. The peak r in the blue doublet member is given 
at the top of each panel. Middle panels show C profiles derived from the doublet solution (equation 3); means and standard deviations 
(sd) measured in the central 150 km s~^ are printed. The actual constant covering factor, C =0.8, is denoted with a dashed line for 
comparison. Right panels show the value of the numerator (N) of the doublet equation divided by C (red), which would be identical with 
the denominator {D, black) for infinite S/N. However, due to the non-linear form of A^, N/C drops below the noise level in D for low r. 
The result is a sytematic underestimation of C, with increasing discrepancy for weaker absorption. 
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